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ABSTRACT 

Several stars orbit around a black hole candidate of mass 3.7 x 10^ Mq, in the region of the 
Galactic Center (GC). Looking for General Relativistic (GR) periastron shifts is limited by the 
existence of a stellar cluster around the black hole that would modify the orbits due to classical 
effects that might mask the GR effect. Only if one knows the cluster parameters (its mass and core 
radius) it is possible to unequivocally deduce the GR effects expected and then test them. In this 
pape r it is shown that the obse rvation of the proper motion of Sgr A*, vsgr A* — (0.4 ± 0.9) km 



s " ( Reid and Bruthaledliool . could help us to constrain the cluster parameters significantly 
and that future measurements of the periastron shifts for at least three stars may adequately 
determine the cluster parameters and the mass of the black hole. 

Subject headings: Gravitation — Galaxy: center — Physics of black holes 



1. Introduction 

GR predicts that orbits about a massive cen- 
tral body suffer periastron shifts yielding rosette 
shapes. However, the classical perturbing effects 
of other objects on inner orbits give an oppo- 
site shift. Since the periastron advance depends 
strongly on the compactness of the central body, 
the detection of such an effect may give informa- 
tion about the nature of the central body itself. 
This would apply for stars orbiting close to the 
GC, where there is a "dark object", the black 



hole hypothesis being the most natural explana- 
tion of the observational data. A cluster of stars 
in the vicinity of the GC (at a distance < 1") has 
been monit ored by ESQ a nd Keck teams for sev- 



eral y ears ( Genzel et al.l ( 2003 al): ISchodel et al 
(l2003l) : iGhez et all (l200l l2004L boOSh l. In 
ticular, Ghez et al. (|2003l ) have reported on ob- 
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servations of several stars orbiting close to the 
GC massive black hole. Among those, the S2 
star, with mass Ms2 — 15 Mq, appears to be a 
main sequence star orbiting the black hole with 
a Keplerian peri od of ~ 15 yrs. This yields 



Ghez et al.l (|2005f ) a mass estimate of Msgr a* — 
3.67 X 10^ Mq within 4.87 x IQ-^ pc, that is the 
S2 semi-major axis. 

Several authors have discussed the possibility of 
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measuring the GR c orrections to Newtonian orbits 
_for Ssr A* (s ee e.g |jaroszvnski"l998', 'l99&', '200C ; 
Fragile & M athewsl[200 0: Rubilar & Eckart 2001 ; 
Weinberg. Miloslavlievic fc Ghea 120051 ). usually 



assuming that the central body is a Schwarzschild 
black hole. However, since black holes generally 
rotate, and there is no reason why they should 
not be rotating fast, the Kerr metric should be 
used instead. Not only stellar mass black holes 
but also supermassive black holes are believed 
to be spinning. Indeed, X-ray observations of 
Seyfer t galaxies, rn i croqu a sars and bina r y sys - 



tems (iFabian et al.l (Il995l): iTanaka et al.l (|l995l ): 



iFabian et all (|2000f i nFabianl (2005( 1 and refer 



ences therein) show that the data could be ex- 
plained by a rotati ng black hole model (see e.g. 



Zakharov &: Repin I (l 2003bl.Q): 'Zak harov eFar 



(|2003ar ) and lZakharov fc Rcpin (2001)). Further, 
supermassive black holes at the center of QSOs, 
AGNs and galaxies show beamed jet emission im- 
plying that they have non zero angular momen- 
tum. Hence, Kerr black holes may be fairly com- 
mon in nature. The relatively short orbital period 
of the S2 star encourages a search for genuine GR 
effects like the orbital periastron shift. Quite pos- 
sibly, more suitable stars, close to the GC black 
hole, will be found in the future. Bini et al. (I2OO5I ) 
studied the GR periastron shift around Sgr A* and 
estimated it for various solutions belonging to the 
Weyl class, including the Schwarzschild and Kerr 
black holes. However, they did not take into ac- 
count the presence of a stellar cluster, which could 
in principle be sizable. 

The purpose of this paper is to try to find lim- 
its for the extent and density of the cluster about 
Sgr A* and if those limits yield a measurement of 
its spin. 

Clearly, a thorough knowledge of the cluster 
mass and density distribution is necessary to be 
able to infer the mass and spin of the black hole 
at the GC by measuring the periastron shift and 
subtracting the Newtonian shift. Unfortunately, 
the star cluster parameters are poorly known. 
However, the measure of the Brownian motion of 
the central black hole due to the surrounding mat- 
ter may be used to constrain the black hole to 



^We remark that the star cluster we are considering around 
the central black hole might contain not only normal stars 
but also white dwarfs, neutron stars and/or stellar mass 
black holes. 



cluster mass ratio0 The latest observations of the 
Sgr A* proper motion, vagr a* = (0.4 ± 0.9) km 



s (jReid and Bruthalerj 120041 ) . is much tighter 



than the earlier one of 2— 20 km s ^ (see Reid, 
Readhead, Vermeulen et al. 

For a test particle orbiting a Schwarzschild 
black hole of mass A/bh, t he pe riastron shift is 
given by (see e.g. Weinberg, Il972| ) 



A0£ 



6ttGMbh 



3{18 + e^)TTG^MlH 
2d2(l -e2)2c4 



(1) 



d and e being the semi-major axis and eccentric- 
ity of the test particle orbit, respectively. For a 
rotating black hole with spin parameter a = |a| = 
J/GMbh, the space-time is described by the Kerr 
metric and, in the most favorable case of equa- 
torial plane motion (a.v = 0), the shift is given 
by (B oyer and Price [l965l . but see also Bini et al. 
20051 for more details) 



AcPk ^ A(j)s + 



(2) 

which reduces to eq. ([!]) for a ^ 0. In the more 
general case, a.v ^ 0, the expected periastron shift 
has to be evaluated numerically. 

The expected periastron shifts (mas/revolution), 
A(j) (as seen from the center) and A(j>E (as seen 
from Earth at the distance i?o — 8 kpc from 
the GC), for the Schwarzschild and the extreme 
Kerr black holes, for the S2 and S16 stars turn 
out to be A^'^^ ^ 6.3329 x 10^ and 6.4410 x 10^ 
and A(/)|;2 = 0.661 and 0.672 respectively, and 
A0^i6 = 1.6428 X 10^ and 1.6881 x 10^ and 
A0|i6 = 3.307 and 3.399 respectively. Recall 
that 

d{l + e) 



Acl)i 



Rr, 



-A(t)s,K 



(3) 



Notice that the differences between the periastron 
shifts for the Schwarzschild and the maximally ro- 
tating Kerr black hole is at most 0.01 mas for the 
S2 star and 0.009 mas for the S16 star. In or- 
der to make these measurements with the required 



^ Other methods for estimating the black hole parameters 
(i.e. mass and angular momentum) based on gravita- 
tional retrolensing have been pr oposed. 



For more de- 



tails o n this topic w e refer to |De Paolis et all ||2003 
[20041) ; IZakharoy et al] ll2005bf) : iDe Paolis et al.l ||2005D 
IZakharov et al. I l l2005al ) and reference therein. 
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accuracy, one needs to know the S2 orbit with a 
precision of at least 10 /J,as. 

There is a proposal to improve the angu- 
lar resolution of VLTI with the PRIMA facil- 
120031); 



model density profile ijBinnev fc TremaindllQSTT) 



itv (IRottgering et al 



Delplancke et al 



(|2003[ ): IQuirrenbach ( 2003f ) but see also the re- 
lated web-sit^, which, by using a phase refer- 
enced imaging technique, will get ~ 10 /ias an- 
gular resolution. Hence, at least in principle, the 
effect of a maximally rotating Kerr black hole on 
the periastron shift of the S2 star can be distin- 
guished from that produced by a Schwarzschild 
black hole with the same mass. 

The plan of the paper is as follows: In the next 
section we briefly discuss the effect of a central 
star cluster on the periastron advance. In Section 
3 we use the Sgr A* Brownian motion to constrain 
the black hole to star cluster mass ratio. Then we 
consider whether the detection of the spin of the 
black hole from the periastron shift of the S2 star 
is possible, once the cluster density and size have 
been adequately constrained. In Section 5 we show 
how future measurements of the periastron shifts 
for at least three stars close to the GC black hole 
may be used to estimate the black hole mass and 
the star cluster mass density distribution. In the 
next section we consider what the observational 
requirements would be for adequate determination 
of the cluster parameters to be able to resolve the 
Kerr effect. Finally, in section 7, we present some 
concluding remarks. 

2. Retrograde shift due to a central stellar 
cluster 

The star cluster surrounding the central black 
hole in the GC could be sizable. At least 17 mem- 
bers have been ob served within 15 mpc up to now 
( Ghez et al. Iliooi). However, the cluster mass and 
density distribution, that is to say its mass and 
core radius, is still unknown. The presence of this 
cluster affects the periastron shift of stars orbiting 
the central black hole. The periastron advance 
depends strongly on the mass density profile and 
especially on the central density and typical length 
scale. 

We model the stellar cluster by a Plummer 



•^http: / /obswww.unigc.ch/PRIMA/homc/introduction. 



PCL{r) = pof{r) , with /(r) 



1+ - 



21 -"/2 



(4) 



where the cluster central density po is given by 

McL 



Po 



J^""- ATTr^f{r) dr 



(5) 



RcL and Mcl being the cluster radius and mass, 
respectively. According to dynamical observations 
towards the GC, we require that the total mass 
M(r) = Mbh + McL{r) contained within r ~ 5 x 
10-3 pc is M ~ 3.67 X 10^ Mq. Useful information 
is provided by the cluster mass fraction, Xcl = 
Mcl/M, and its complement, Xbh ~ ^ ^ Xcl- 
As one can see, the requirement given in eq. ([5|) 
implies that M{r) Mbh for r ^ 0. The total 
mass density profile p{r) is given by 



p(r) = AsHM,5(3)(r)+po/(r) 
and the mass contained within r is 



(6) 



M{r) = XbhM+ / 47rr'2po/(/) dr' . (7) 



In Figure [T] we show the cluster mass density pro- 
file pcL{r) as given by eq. (j4|), for selected values 
of Xbh- The total mass M{r) enclosed within the 
radius r is also shown in Figure [2] In both Fig- 
ures, solid, dotted and dashed lines correspond to 
Xbh ~ 0, 0.7, 0.9, and we have assumed Tc = 3 
mpc (thick lines) and Vc = 5.8 mpc (thin lines). 

The Newtonian gravitational potential $Ar at a 
distance r due to the mass contained within it can 
be evaluated as 

M(r') 



$jv(r) = -G f 

J r 



dr' 



(8) 



In FigureEl the gravitational potential <I>jv(r) due 
to the mass density distribution in eq. ^ is given 
for selected values of Xbh- 

According to GR, the motion of a test parti- 
cle can be fully described by solving the geodesic 
equations. Under the assumption that the matter 
distribution is static and pressureless, the equa- 
tion of motio n of t he test particle becomes (see 
e.g. Weinberg Il972h 



dv 
dt 



-V($ 



N 



4v(v • V)$7V - U^V$Ar . 

(9) 
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Fig. 1. — The cluster mass density profile is shown 
for different values of \bh- Solid, dotted and 
dashed hnes correspond to \bh = 0, 0.7, 0.9, 
respectively. Thick, red lines have been obtained 
for Tc = 3 mpc with the same values A bh as given 
above. 




For a spherically symmetric mass distribution 
with a density profile given by eq. ([4]) and for a 
gravitational potential given by eq. ([8|), the pre- 
vio us rel ation becomes (see for details Rubilar et 



Fig. 2. — The mass enclosed within the distance 
r is shown for different fractions \bh of the total 
mass M which is contained in the central black 
hole. Solid, dotted and dashed lines correspond 
to Xbh = 0, Xbh = 0.7 and Xbh = 0.9, re- 
spectively. Note that the case corresponding to 
Xbh = is not real i stic a s shown by some obser- 
vations (|Shen et al.l (|2005l) ). 



al. I2OOII) 



dv 

lit 



GM{r) 



4$ 



N 



4v(v ■ 



(10) 

r and v being the radius vector of the test particle 
(with respect to the center of the stellar cluster) 
and the velocity vector, respectively. Once the ini- 
tial conditions for distance and velocity are given, 
the orbit of a test particle can be found by solv- 
ing the set of ordinary differential equations in eq. 
numerically. 
Now consider the S2 star, which is moving 
around the central distribution of matter on an 
elliptic orbit of semi-major axis d and eccentric- 
ity e in the Newtonian approximation. We take a 
frame with the origin in the GC, X-Y plane on the 
orbital plane and the X axis pointing toward the 
periastron of the orbit. Hence, we can choose the 
Newtonian initial conditions to be (see e.g. Smart 



= d(l + e) 
= , 



(11) 



and 



lGM{rl) 



(12) 



For the S2 star, d and e given in the literature 
are 919 AU and 0.87 respectively. They yield the 
orbits of the S2 star for different values of the black 
hole mass fraction A bh shown in Figure ID The 
Plummer model parameters are a = 5, core radius 
Tc — 5.8 mpc. Note that in the case of Xbh = 1, 
the expected (prograde) periastron shift is that 
given by eq. H]), while the presence of the stellar 



''We would like to mention that the dynamical state of the 
region around Sgr A* is known to be complex, with a sig- 
nificant population of young stars of unclear origin making 
the assumption of an undisturbed spherical cluster likely 
uncorrect. Considering the effects caused by a non spheri- 
cally symmetric mass distribution makes the passage to an 
equation similar to eq. mOI I not analytically solvable. This 
problem will be addressed in a subsequent work. 
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cluster leads to a retrograde periastron shift. For 
comparison, the expected periastron shift for the 
S16 star is given in Figure O In the latter case, 
the binary system or bital parameters were taken 
from Schodel et al. (|2003l) assuming also for the 
S16 mass a conservative value of ~ 10 M0. 

In Figure [S] the S2 orbital shift is given as 
a function of the stellar cluster core radius r^, for 
different power law index values (a = 4 dashed 
line, a — 5 dotted line and a = 6 solid line). 
In the left panel, the black hole mass fraction is 
Afl ff — .8 in order to compare with Rubilar et 
al. (|200ll ) results, while the right panel shows the 
^BH = 0.99 case. Note that for extremely com- 
pact clusters, A$ is quite small. The same is true 
for large enough core radii, corresponding to mat- 
ter density profiles almost constant within the S2 
orbit. 

Figures [4] and [6] show that the expected S2 peri- 
astron shift depends strongly on the total mass of 
the cluster. In particular, the shift due to the clus- 
ter is opposite in sign (retrograde motion) to the 
relativistic effect due to the black hole in the GC. 
Moreover, for each value of the cluster mass and 
power law index, there exist two density profiles 
(corresponding to two particular core radii) which 
have total shift almost zero, implying that the pe- 
riastron advance due to the cluster is equal (in 
magnitude) to the periastron shift due to the black 
hole. A numerical analysis shows that the transi- 
tion from a prograde shift (due to the black hole) 
to retrograde shift (due to the extended mass) 
occurs at Xbh ^ 0.9976, 0.9986 and 0.9990 for 
a = 4, 5 and 6, respectively. This means that 
a small fraction of mass in the cluster drastically 
changes the overall shift. 

We would like to note that the assumption of 
the Plummer model to describe the mass density 
distribution of the stellar cluster around the cen- 
tral black hole is an oversimplification. Indeed, 
one expects that in presence of a central black 
hole, the stellar profile should follow a Bachall- 
Wolf law with density distribution pdr) oc r~'^/'^ 
(|Bahcall fc Woljll977t iBinnev &: Tremainel IT987h 



at least up to th rn, where vh = GMBH/f^i — 
0.5 pc is the radius of the black hole influence 
sphere. In the following, we call tr the distance 
(<C rn) up to which the cluster mass density fol- 
lows the Bachall-Wolf law. 

In order to study the effect of such a cusp on the 



expected S2 periastron shift, we consider three dif- 
ferent cases a) the cusp is entirely contained within 
the S2 periastron distance Rs2 (i-e. fn < -^52), 
b) the cusp extends beyond the S2 periastron dis- 
tance (thus making the S2 star move in a mass 
gradient) and c) the stellar density profile follows 
a cusp law up to the distance fn from the cen- 
ter and a Plummer law for r > fn- In cases a) 
and b) all stars are in a cusp density profile. In 
any case we require that the total mass enclosed 
within 4.87 x lO'^ pc is M ~ 3.67 x 10^ Mq. 

In case a), the total S2 periastron shift is just 
the sum of the shift due to the black hole and the 
shift caused by the stellar cusp (that contributes 
with the same sign). Hence, the S2 shift turns out 
to be A$ ~ 0.17 degree per revolution. 

In case b ), by requiring that the total mass en- 
closed within 4.87 X 10-3 pc is M ~ 3.67 x 10^ M©, 
we find that the dependence of the cusp mass and 
the induced S2 periastron shift on fn vanishes. In- 
deed, in Figure[7l we give the mass enclosed within 
the distance r for different values of Xbh- Sohd, 
dotted and dashed lines correspond to Xbh = 0.1, 
Xbh = 0.5 and Xbh = 0.9, respectively. Figure[8] 
shows the expected S2 periastron shift as a func- 
tion of Xbh- As noted before, the shift due to the 
cluster is opposite in sign with respect to that due 
to the black hole. Moreover, for Xbh — 0.998 the 
total shift turns out to be zero since the contribu- 
tions of the black hole and the cluster cancel out. 
It is noticing that, since in the case of cusp profiles 
the density gradient is larger than in the case of a 
usual (a = 5) Plummer model, the value of the S2 
periastron shift gets generally larger values. Only 
if the Plummer core radius is around 2 mpc the 
resulting S2 periastron shifts are comparable in 
both cases. We have then considered the super- 
position of a Plummer model and a Bahcall-Wolf 
profile (case c) extended up to fn such as the cusp 
density at fu equals that of the Plummer model 
at the same distance. Here, if th ^ Rs2 the S2 
periastron shift will be practically equal to that 
caused by the Plummer model (see right panel in 
Fig. 6) since in this case the cusp will have a minor 
influence. On the contrary, for an extended cusp 
[fn S> Rs2)^ the cusp effect on the S2 periastron 
shift will dominate reconciling with case b. 

As a last point, we mention that we have also 
considered the effect due to an extrapolation of 
the observed stellar density profile - the inner- 
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Fig. 4. — Post Newtonian orbits for different values of the black hole mass fraction \bh are shown for the S2 
star (upper panels). Here, we have assumed that the Galactic central black hole is surrounded by a stellar 
cluster whose density profile follows a Plummer model with a = 5 and a core radius ~ 5.8 mpc. The 
periastron shift values in each panel is given in arcseconds. 




Fig. 5. — The same as in Figure [4] but for the S 16-Sgr A* binary system. In this case, the binary system 
orbital parameters were taken from Ghez et al. ( 20051 ) assuming for the S16 mass a conservative value of 
~ 10 Mq. 
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Fig. 6. — The expected S2 periastron shift as a function of the stehar cluster core radius is shown. Here we 
have assumed a Plummer density profile for the stellar cluster. Dashed, solid and dotted lines correspond 
to a = 4, 5 and 6, respectively. The black hole mass fraction has been fixed to Xbh = 0.8 (left panel) and 
\bh = 0.99 (right panel), respectively. Note the existence of a maximum approximately corresponding to 
the S2 semi-major axis. 
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Fig. 3. — The gravitational potential at distance 
r as due to the mass M{r) is shown for different 
fractions Xbh of the total mass M . Solid, dotted 
and dashed lines correspond to Xbh — 0, Xbh = 
0.7 and Xbh = 0.9, respectively. Thick red lines 
have been obtained for Vc — 5 mpc while thin black 
lines are for Vc = 5.8 mpc. 



most point of which is the S2 star at a di s tance o f 
0.1"- within i?s2- Following'GenzeietalJ (l2003_b|) 
and assuming a cusp stellar density profile, we 
find that the enclosed mass is in the range 30-300 
Mq (for a constant mass density or a power law 
with index 7 = 1.4). Therefore, the cusp effect on 
the S2 periastron shift is negligible since the cor- 
responding Xbh is always greater than 0.99992. 
However, we caution that the case under inves- 
tigation in__the_gresent_2a£er is different with re- 
spect to Genzel et al. ( 2003 bf ) since we are assum- 
ing that a fraction of the mass contained within 
Rs2 may be in a stellar cluster. Hence, the clus- 
ter mass content may be larger, thus providing a 
stronger effect on the S2 periastron shift. 

3. Tightening mass limits of Sgr A* 



We know that the mass of Sgr A* within the S2 
orbit is 3.67 x 10® M© to a high accuracy. Though 
there is nothing definite known about the mass 
distribution, there is strong reason to believe that 
there is a black hole of several solar masses, possi- 
bly surrounded by a significant cluster. In princi- 
ple the cluster mass could dominate over the black 
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Fig. 7. — The mass enclosed within the distance 
r is shown for different fractions Xbh of the to- 
tal mass M contained within the S2 orbit. Solid, 
dotted and dashed lines correspond to = 0.1, 
Xbh = 0.5 and Xbh = 0.9, respectively. The stel- 
lar cluster is assumed to follow an r~^/^ density 
profile. 

hole, be comparable to it or be dominated by it. 
That there is a cluster is highly likely on account 
of the large number of stars observed near Sgr A* . 
Though these lie outside the S2 orbit, many stars 
so far unseen probably do lie within the orbit as 
well. In this section we use current data on the 
Brownian motion of Sgr A* and the evaporation 
time for the putative cluster to put limits on the 
cluster mass and hence on the black hole mass. 

Chatterjee, Hernquist and Loeb ( 2002t ) have de- 
veloped a simple model to describe the dynamics 
of a massive black hole surrounded by a dense stel- 
lar cluster. The total force acting on the black hole 
is separated into two independent parts, one of 
which is the slowly varying force due to the stellar 
ensemble and the other the rapid stochastic force 
due to close stellar encounters. In the case of a 
stellar system with a Plummer distribution, the 
motion of the black hole is similar to that of a 
Brownian particle in a harmonic potential. Thus 
the black hole one-dimensional mean-square veloc- 
ity is given by 



n 2 GAIcLm^ 

< V >= 

9 rcMBH 



(13) 



Fig. 8. — The expected S2 periastron shift as a 
function of the mass ratio parameter Xbh- Solid 
and dashed lines correspond to the S2 shift due to 
the black hole and to the stellar cusp, respectively. 
We note that the shift due to the stellar cusp is 
independent on the fn value, that, in this case, 
has been assumed to be larger than the S2 semi- 
major axis (case b). 

where it has been assumed that the cluster is com- 
posed of objects with equal mass, m*. For a Plum- 
mer {a — 5) stellar cluster, the total mass within 
R is 



M{R) = M_ 



BH 



{R'^+rlfl^ 



(14) 



Since < > is less than a certain maximum 
value < v1 >max, from eqs. (fT^ and pi)) one 
obtains 



Mbh > M{R) |l 



G{R^ 



. . 

the right hand side corresponding to a minimum 
black hole mass, as constrained by the Brownian 
motion of the central black hole. In Figure [9] the 
minimum black hole mass allowed by the Brown- 
ian motion of Sgr A* is given as a function of the 
stellar cluster core radius, for two different proper 
motion velocities of the black hole: 1.3 km s~^ 
(dashed hues), and 2 km s~^ (dotted hues). The 
total mass contained within R = 0.01 pc of Sgr A* 
has been taken to be M ~ 3.67 x 10^ Mq. 

Chatterjee, Hernquist and Loeb ( 20021 ) derived 
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Fig. 9. — The minimum black hole mass allowed 
by the Brownian motion of Sgr A* is given as a 
function of the stellar cluster core radius for the 
different black hole proper motion velocities. We 
assume that a total mass M ~ 3.67 x 10^ is 
contained within Rs2 = 4.87 mpc. Dashed and 
dotted lines have been obtained for velocities of 
1.3 km s""'^ and 2 km s~^, respectively. For each 
given curve only the region above it is allowed. 

an evaporation time for a cluster, but concen- 
trated on the large scale cluster Tc — 10 pc about 
Sgr A*, and hence assumed that Mcl ^ Mb//. 
On the other hand Rauch and Tremaine ( |1996[ ) 
and Mouawad et al. ( 2005f l consider only the re- 
gion interior to the orbit of S2 and assume M-cl 
< Mbh- 

We need to allow for all possibilities while con- 
sidering the cluster interior to the orbit of S2, in- 
cluding Mcl — y^BH- For this purpose we con- 
sider a cluster of core radius Tc and mass Mcl = 
M — M BH ■ We now need to obtain the generaliza- 
tion o f the formula of Chatterjee, Hernquist and 
Loeb (|2002l ) for the median relaxation time in this 
more general situation. For this purpose, as usual, 
we assume that the cluster consists of components 
of the same mass and evaluate the crossing 
time in the usual way to obtain the general me- 
dian relaxation time 



Tr = 



0.14(1.3 TcM) 



3/2 



It is easy to verify that in the approximation Mcl 
, » M-BH we recover the fo rmula of Chatterjee, 
Hernquist and Loeb (|2002r ) and in the approxi- 
mation Mcl ^ M-B H recover the formula of 
Rauch and Tremaine ( 19961). The evaporation time 



is then T^^ap ^ 300 % (jBinnev fc Tremainelll987 



GMcLinJog{QAM /m^ 



(16) 



p. 525). 

One can assume different "reasonable" values 
of the time that the cluster would have been in 
existence and hence use the evaporation time to 
further limit the black hole mass in the GC. It 
is clear that 10^ years = 0.1 Gyr is less than the 
minimum value that could be regarded as reason- 
able, 1 Gyr is more reasonable and 10 Gyr is likely 
^ to be a good value to take. The results are given 
i°m Table [U for m, = 1 M©. Note that the tight- 
est bound gives a very stringent upper limit of 
9 X IO^Mq on the cluster mass. Also note that 
the value decreases if the average is taken to 
be larger. 

4. The spin of the black hole 

The periastron shift is the net contribution of 
the relativistic retrograde shift due to the black 
hole and the Newtonian prograde shift due to 
the surrounding cluster. Obviously, if the pe- 
riastron advance due to the stellar cluster were 
known, the contribution of periastron advance due 
to the black hole could be obtained by subtracting 
from the measured quantity. The question arises 
whether the information obtained would be ade- 
quate to obtain both the black hole mass and spin 
parameters. Though we can put reasonably sharp 
bounds on the stellar cluster about the black hole, 
is it good good enough for our purpose? If so, we 
could use eq.([2]) to obtain the spin of the black 
hole for different values in the possible range for 
the periastron shift. It is easy to see from Fig[6] 
that for \bh — 0.99 and allowing for the max- 
imum range of unknown values of a and rc the 
1.8 X 10^ < -A0 < 4.7 X 10^ or 1.9 x 10"^ < 
—A(j)E < 4.7 X 10^^. For the sharpest limit ob- 
tained, a = 5, and Brownian motion 1.3 km s~^, 
Xbh = 0.975, we find that A<j) ~ -4.47 x 10^ 
or A(pE — —4.7 X 10^^. For Brownian motion 
2.0 km s-\ Xbh = 0.964, A0 ~ -5.75 x 10^ or 
A4>E — —6.0 X 10^'^. This is a factor of 5 less than 
the effect of the spin. Hence this method cannot 
be used to determine the spin. For this we need 
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Tovap(Gyr) 


rj-3(mpc) 


\1.3 


r^"(mpc) 


\2.0 


0.1 


0.87 


0.762 


1.28 


0.645 


1 


2.11 


0.919 


2.61 


0.876 


10 


4.12 


0.975 


5.27 


0.964 



Table 1: The cluster core radius and minimum black hole mass fraction Xbh for the limits obtained by 
< >max= 1-3 and 2.0, for T^vap = 0.1, 1 and 10 Gyr. 



the cluster parameter values, rather than upper 
limits for them. Alternatively, one would need to 
rely on the retrolensin g method suggested earlier 



(|De Paohs et al.ll2005l : ?) 



5. Determination of cluster parameters 

Using the stronger (Icr) limit of 1.3 km and 
the weaker {2a) limit of 2.0 km s^^ to limit the 
Brownian motion of Sgr A* for our calculations 
and evaporation times of 1 and 10 Gyr for the clus- 
ter, we obtained the minimum black hole mass. 
For the stronger limits on the Brownian motion 
and the evaporation time, it is 3.579 x 10^ M0 
corresponding to a Xbh — 0.975 for q = 5. Our 
numerical analysis shows that the transition from 
a prograde shift (due to the black hole) to a ret- 
rograde shift (due to the extended mass assumed 
to be distributed with a Plummer density profile) 
occurs at Xbh ^ 0.9976, 0.9986 and 0.9990 for 
a = 4, 5 and 6, respectively. Hence, even a small 
cluster around the central massive black hole lim- 
its the possibility to observe and use the periastron 
shift of the S2 star. 

Since we have modeled the star cluster density 
profile by a Plummer model, the periastron shift 
contribution due to the stellar cluster depends on 
three parameters: the central density po (or equiv- 
alently Xbh)', the core radius rc', and the power- 
law index a. This degeneracy in the determination 
of the stellar cluster parameters is due to the mea- 
surement of the periastron shift of a single star. 
This is easily seen by inspecting Figure 1101 which 
has been obtained for illustrative purposes for the 
S2 star by setting Xbh = 0.99 and varying both 
the core radius Vc and power-law index a for the 
star cluster density profile. Each contour line cor- 
responds to a given S2 periastron shift in units of 
degrees. To solve the parameter degeneracy and 
determine the stellar cluster parameters (by study- 



ing the periastron advance effect), the periastron 
shifts for at least three different stars have to be 
measured with sufficient accuracy. Consider, for 
example, the S16 star having an orbital period of 
~ 36 yr and eccentricity e ~ 0.97. Measuring its 
periastron shift and comparing with the S2 result 
will give much tighter information about the stel- 
lar cluster parameters. From Figure [11] it is evi- 
dent that there are regions (intersections between 
dashed and solid lines) in the a-rc plane for which 
one measures values of the periastron shift for the 
S2 and S16 stars. Obviously, there could be (as 
yet unobserved) stars with orbit apocenters com- 
parable to S2, but with different eccentricities (for 
example larger than 0.87) or stars closer to the GC 
black hole than S2 or S16 stars. Monitoring their 
orbits and measuring their periastron shifts will 
be extremely helpful in reconstructing the cluster 
density profile. As an example, in Figure [12] we 
compare the expected S2 periastron shift (solid 
lines obtained for Xbh = 0.99) with the perias- 
tron shift of a star whose orbit has an eccentricity 
of ~ 0.87 and semi-major axis 3 times smaller than 
that of S2. 

As is evident from Figures 8-10, one can ob- 
tain estimates of the rc, ol and Xbh, provided that 
three stars have been observed to sufficient accu- 
racy. Assume that we have adequate accuracy of 
observation to see periastron shifts of 10~^'^ mas, 
which is the value required to see the relativistic 
periastron shift. To what accuracy have we lim- 
ited the cluster parameters? To determine this, we 
could just vary Xbh for a given Vc- The effect of 
this change would be less than the effect of chang- 
ing Tc and Xbh- As such, if we want to know how 
accurately the cluster parameters are determined, 
we need to calculate the maximum change in Vc 
along with the change in Xbh, as allowed by the 
Brownian motion limit. By varying Xbh by 10~^ 
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Fig. 10. — The expected S2 periastron shift for 
As/f = 0.99 for different values of both the core 
radius Tc and power-law index a for the consid- 
ered Plummer density profile. Each contour line 
corresponds to a given S2 periastron shift in degree 
units. Note that a degeneracy occurs since there 
exist different values of the power law index and 
core radius corresponding to the same periastron 
shift. 

and Tc maximally we find that we get the required 
accuracy, for evaporation rates from 1 to 10 Gyr 
and Brownian motion of 1.3 to 2.0 km . With 
this accuracy we should also be able to separate 
out the classical periastron shift of the stellar clus- 
ter and the relativistic effect of a maximal (and 
even slightly less than maximal) Kerr black hole. 
With better accuracy we should be able to get an 
estimate, or at least an upper bound, for the black 
hole spin as well. The question now is, what is re- 
quired to achieve this accuracy of observation of 
the periastron shift of three stars? This point is 
discussed in the next section. 

6. Observational requirements for deter- 
mination of black hole spin and cluster 
parameters 

In the near future, observations using large di- 
ameter telescopes in combination with adaptive 
optics may allow us to reach the angular reso- 
lution needed to measure the periastron shift of 



Fig. 11. — The same as in Figure 8. Dotted lines 
show contours for the S16 star. If the periastron 
shifts of both stars will be measured in the future 
the intersection between the corresponding con- 
tour lines will give information about the central 
stellar density profile. 



stars close to the GC back hole. Consider for ex- 
ample an instrument with an angular resolution 
Ac&yi and assume that the relative position of stars 
can be determined to about 1/e of the achieved 
angular resolution, i.e. the position accuracy is 
A$p ~ A<l>A/e. The positional accuracy can be 
increased by a factor V^, if N reference stars are 
used. In this case, the maximum positional acc u- 
racy is simply given by (|Rubilar fc Eckartj|200l[ ) 



A$, 



(17) 



It follows that if the periastron position of a 
star shifts by an amount A^^; (as observed from 
Earth) , to obtain the desired accuracy we need at 
least that A^p ~ A^^. In this case, the mini- 
mum number of reference stars can be determined 
once both the instrument angular and positional 
accuracies are known, i.e. 



^ ni.ax — 



/A$j 



V A$i 



(18) 



As an example, the LBT interferometer has 
angular resolution A$^ ~ 30 mas and the rela- 
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Fig. 12. — The same as in Figure 8. Dotted lines 
show contours for a star with orbit of ~ 0.87 (the 
same as the S2 star) but semi-major axis 3 times 
smaller with respect to the S2 one. If the peri- 
astron shifts of both stars will be observed in the 
future the intersection between the corresponding 
contour lines gives information about the central 
stellar density profile. 



five position of stars is conser vatively estimated 
to be about 1/30 of that value ( Rubilar k, Eckarj 



20011 ). Therefore, to measure the periastron shift 



with adequate accuracy to see the least shift for 
the cluster parameters allowed, and thereby detect 
the relativistic shift, we need to be ~ (6.6 x 
10~^/2 X 10~^)^ or 10^ reference stars, which au- 
tomatically provides the accuracy required to see 
the maximal Kerr (spin) effect. For PRIMA (see 
^ttp://www. eso.org/projects/vlti/instru/prima/index. 
the relative positional accuracy is planned to be 
~ 10 /ias. As such, we only need a single reference 
star. 

7. Concluding remarks 

We have used the fact that the stellar cluster 
close to the central black hole seems spherically 
symmetric to limit the Brownian motion of Sgr 
A* to be the observed proper motion. We have 
taken the stronger (Icr) limit of 1.3 km s~^ and 
the weaker {2a) limit of 2.0 km s~^ for our cal- 
culations. We also used evaporation times of 1 



to 10 Gyr for the cluster, appropriately modified 
to incorporate the gravitational well due to the 
E black hole, to put further constraints on the clus- 
ter mass. The results of our calculations show 
E that the stellar periastron shifts due to the clus- 
E ter, even limited to the extent considered, may 
E totally swamp not only the Kerr (spin) effect but 
E also the Schwarzschild effect. However, the discus- 
E sion focused on the observations for a single star, 
E S2. By modelling the star cluster density profile 
E with a Plummer low, the periastron shift contri- 
E bution due to the stellar cluster depends on three 
E parameters: the central density pq (or equivalently 
E Xbh)i the core radius r^, and the power-law index 
E a. Consequently, with observations of three stars 
I'Hve should be able to determine the cluster parame- 
ters adequately. We have addressed the question 
of what is required to obtain the desired accuracy 
for observing the relativistic effect. It turns out 
that we need about 10^ reference stars with the 
LBT interferometer. With the accuracy expected 
of PRIMA, it should be enough to use only one 
reference star. 
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